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C^. in the QCD plasma, we obtain effective amplitudes for the ultrasoft colour fields 

^^ . carrying momenta of order g^T. These amplitudes are of the same order in g as 

the hard thermal loops (HTL), which they generalize by including the effects of 

^ . the collisions among the hard particles. The ultrasoft amplitudes share many of 

» . the remarkable properties of the HTL's: they are gauge invariant, obey simple 

Ward identities, and, in the static limit, reduce to the usual Debye mass for the 

electric fields. However, unlike the HTL's, which correspond effectively to one-loop 
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bare perturbation theory. By solving the linearized Boltzmann equation, we obtain 

a formula for the colour conductivity which accounts for the contributions of the 

hard and soft modes beyond the leading logarithmic approximation. 
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1 Introduction 

In recent years, kinetic theory has proven to be a powerful tool to construct effective 
theories for the soft fields in ultrarelativistic plasmas. Thus, the effective theory at the 
scale gT follows from a coUisionless kinetic equation, of the Vlasov type |I[]. The effective 
theory at the scale g^T is generated by a Boltzmann equation which includes a collision 
term for colour relaxation [^ ||, §, ^. (Here, T is the temperature, and g is the cou- 
pling constant, assumed to be small.) The kinetic description relies on the separation 
of scales between single-particle and collective excitations. This allows for kinematical 
approximations which, like the relevant scales themselves, are controlled by powers of 
g. By using these approximations, kinetic equations have been rigorously constructed 
from the quantum equations of motion [|l], ||, ^, |^, thus providing justification for nu- 
merous previous works using ad hoc transport equations inspired by classical physics 
10, |, ^ [10], [11], |12|, H, |T3|, [I^. Previous attempts to derive these equations |]T5| generally 



failed to recognize the proper separation of scales which turns out to be essential in order 
to control the gauge invariance of the approximations involved. 

The single-particle excitations of the QCD plasma are hard transverse gluon^ with 
typical momenta k r^ T. The so/t fields are colour fields A'^{x) with momenta of the order 
gT or less. When acting on the hard particles, these soft fields induce longwavelength 
(A ^ ^/gT) collective excitations, with A much larger than the mean interparticle distance 
f ~ 1/T |l^, |l^, 1^. In the framework of the kinetic theory, these excitations are described 
by a colour density matrix 6Nab(k;X) to which the soft fields ^^(a;) couple via kinetic 
equations. By solving these equations, one can express SNat as a functional of the fields 
A'^. The corresponding colour current (with f'^ = (l,k//i;)) : 



d^A; 



^^-^^2gJ-—-v^Ti(T'^6N{k,x)), (1.1) 



(27r) 
acts as a generating functional for the (equilibrium) amplitudes of the soft fields ||1], 0: 



r, = W'Al + lT';^:A^,A^^+... (1.2) 



Here, H^^ = 5"*H^j, is the soft polarization tensor, and the other terms represent vertex 
corrections. These are the amplitudes which define the effective theory for the soft fields. 
When applied to colour excitations at the scale gT f!], ^ |16| , this strategy provides the so- 
called "hard thermal loops" (HTL) |]TB|, |T^, |T7|. It is our purpose in this paper to generalize 
this strategy to the ultrasoft scale g'^T, and construct the corresponding amplitudes. 

Remarkably, to the order of interest the density matrix can be parametrized as: 

6Nab{k,x) = -gWat{x,^r){dNo/dk), (1.3) 



^We consider here a purely Yang-Mills plasma, with no quarks. 



where A^o(^) = '^/{e^^ — 1) is the Bose-Einstein thermal distribution, and W{x,\) = 
Wa{x,v)T"' is a colour matrix in the adjoint representation which depends upon the 
velocity v = h/k (a unit vector), but not upon the magnitude k = |k| of the momentum. 
Then, the kinetic equations are written as equations for Wa{x,\). 

Let us briefly recall the situation at the scale gT. The relevant kinetic equation is 
a non-Abelian generalization of the Vlasov equation [|l|: 

{vD,r'Wt{x,v) = vE'^(x). (1.4) 

It differs from the corresponding Abelian equation, namely (with W{x, v) a fluctuation 
in the electric charge density) 

{v-d^)W{x,^^) = v-E(a;), (1.5) 

merely by the repacement of the ordinary (soft) derivative dx ~ gT by the covariant one 
Dx = dx + igA. Accordingly, the soft gluon polarization tensor derived from eq. (|1.4| ) is 
formally identical to the photon polarization tensor obtained from eq. ( p..5[ ). In addition, 
eq. ( p..4|) also generates, through the covariant derivative, an infinite series of gluon ver- 
tices. These are the HTL's alluded to before, corresponding to one-loop diagrams with 
soft external lines and hard internal momenta [|1^, [1^. Note that the kinetic equation 
( |1.5|) isolates directly the dominant contributions of such diagrams, in a gauge invariant 
manner. 

This close similitude between the response of Abelian and non-Abelian plasmas to 
longwavelength perturbations disappears, however, when going to very soft perturbations, 
where collisions start to play a role. The effects of the collisions depend upon the specific 
excitations one is looking at. To give a crude estimate of these effects, one may use the 
relaxation time approximation, where the kinetic equation is written as 

Wix v) 
{vDxtW,{x,Y) = vE-(x) ^^^, (1.6) 

Tcol 

and Tcol is the typical relaxation time for small off-equilibrium colour fluctuations. Like the 
damping rate 7 for hard quasiparticles |^ |21|, ^, ^, to which it is intimately related 



(see below), the relaxation of colour is dominated by the singular forward scattering 
(i.e., by soft momentum transfers in the collision in Fig. |I]), which yields Tcoi ~ I/7 ~ 
1 / {g'^T \n.{l / g)) [§, |^. Then, eq. ( |1.6|) shows that the effect of the collisions become 
a leading order effect for inhomogeneities at the scale dx ~ g^T, or less. This should 
be contrasted with the case of colourless fluctuations^, for instance fluctuations in the 



'^Note also that, for colourless fluctuatfons, the analogue of eq. (1.6) will generally involve a momentum- 
dependent relaxation time |l2|; see Sec. 4.2 below. 



momentum or the electric charge distributions, where the typical relaxation time is much 
larger, Tei ~ l/{g'^T\iii{l/g)), as it requires large angle scattering |0, |lT], |12 . 

For colour fluctuations at the scale g'^T, it is further convenient to constrain the 
amplitudes of the associated mean fields such as 1^4^ | ~ gT; then the two terms of the 
ultrasoft covariant derivative are of the same order in g (namely dx ~ gA ~ g^T) and, in 
the derivation of the kinetic equations, one can consistently preserve gauge symmetry with 
respect to the background field 0]. There is another reason which makes this constraint 
interesting: |yl^| ~ gT is the typical amplitude of the thermal fluctuations at the scale 
g^T [^. These fluctuations have relatively large amplitudes because of Bose-Einstein 



enhancement, and their dynamics is fully non-linear; as a result, perturbation theory 
breaks down at the scale g^T [|l^, 25|. Moreover, these large amplitude fluctuations make 



it impossible to give a gauge independent meaning to inhomogeneities on scales much 
larger than 1/g'^T. A convenient strategy to deal with this situation is to observe that the 
soft modes can be treated as classical fields, precisely because of their large occupation 



numbers [^ ^ ^, (and references therein). Then, the non-perturbative dynamics can 



be studied via classical lattice simulations of the effective theory for soft fields p8|, |2^, [3 



In order to explicitly construct this theory, however, one needs to go beyond the 
relaxation time approximation (|1.6|) . In fact, eq. ( |1.6|) is inconsistent with gauge symme- 
try, as it leads to a colour current which is not conserved. The correct kinetic equation, 
as derived in [^ Q (see also Refs. 0, ||, |13|)! involves a more complicated collision term, 
which is local in x, but non-local in v (see eq. ( ^.15[ ) below). Since this collision term 
is saturated by soft momentum transfers (it is logarithmically sensitive to all momenta 
q ^ gT), it is useful to isolate the ultrasoft (g ~ g'^T) background fields from the soft 
{g^T ^ q ^ gT) gluons exchanged in the collisions by introducing an intermediate scale fi 
such as g^T <^ /i ^ gT (e.g., fi ~ f7^Tln(l/f7)). Then, the Boltzmann equation generates 
an effective theory for the ultrasoft (g < /i) fields, corresponding to "integrating out" the 
hard and soft (g > fi) fields to leading order in perturbation theory p|, ^ (see also Sees. 
2 and 3.2 below for a discussion of the relevant approximations). The scale /i acts as an 
infrared (IR) cutoff for the collision integral, and as an ultraviolet (UV) cutoff for the 
effective theory, and it must cancel in any complete calculation of ultrasoft correlation 
functions (a cancellation referred to as matching). 

It is our purpose in this paper to study the contribution of hard and soft fields to 
the amphtudes with ultrasoft external fields {ultrasoft amplitudes in brief) by an analysis 
of the solution to the Boltzmann equation. 

In previous applications of the latter — namely, to the calculation of the (transverse) 
colour conductivity to leading logarithmic accuracy [§, ||, |^, |] — , the non-local piece of 
the collision term turned out not to be important. But this was specific to that particular 



approximation, which has ignored all the non-local and non-linear effects in the problem 
(essentially because the drift term v ■ D^ ^ g^T has been neglected as compared to 
7 ~ g'^T\n.{l/ g)] see Sec. 3.4 for more details). In that situation, the ultrasoft amplitudes 
in eq. (|l]2|) collapsed to a single, local quantity, namely the colour conductivity. 

Our intention here is to go beyond this leading logarithmic approximation and 
study the generic ultrasoft amplitudes generated by the Boltzmann equation for colour 
relaxation. This includes, in particular, the non-local effects in space and time, as governed 
by the drift term v ■ D^, and also the non-local effects in v coming from the collision term. 
Because of the latter, the Boltzmann equation cannot be exactly solved in general. Thus, 
we will not be able to provide full expressions for the ultrasoft amphtudes except in some 
simple limits (cf. Sec. 4.2). Still, many of the important properties of these amplitudes 
can be inferred from an analysis of the Boltzmann equation (cf. Sees. 3.1, 3.3 and 4 
below). Moreover, some formal solutions can also be obtained, by iterations, and this is 
specially useful for comparison with diagrammatic perturbation theory (cf. Sees. 3.2 and 
3.3). The Boltzmann equation accomplishes a non-trivial resummation of the perturbative 
expansion, made evident by the derivation of this equation from quantum field theory |^, 
where the collision term has a direct diagrammatic interpretation (see Sec. 2 below for a 
short review of this derivation). 

Let us briefly enumerate here the main properties of the ultrasoft amplitudes, to be 
derived below in this paper: For generic external momenta of order g'^T, these amplitudes 
are of the same order in g as the HTL's, which they generalize by taking into account the 
effects of the collisions. The ultrasoft amplitudes share many of the remarkable properties 
of the HTL's: i) they are non-perturbative, in the sense that in the kinematical regime of 
interest (tu <^ p ^ g^T), they are as large as the corresponding tree-level amplitudes; ii) 
they are gauge-fixing independent, as they involve only collisions among on-shell, hard, 
transverse gluons, in) they satisfy simple Ward identities, which express the conservation 
of the colour current; iv) in the static limit cj — > 0, they reduce to the usual Debye mass 
term m|, = g'^NT'^/3 for the electric fields. That is, to the order of interest, collisions do 



not modify Debye screening (see also Ref. |31], and Sec. 3.3 below for more details). 



On the other hand, there are also significant differences with respect to the HTL's: 
i) The ultrasoft amplitudes have no Abelian counterpart: in QED, the effects of the 
collisions become important only at the scale e^T (see also Sec. 2 below), ii) Unlike the 
HTL's, which correspond to one-loop Feynman graphs, the ultrasoft amplitudes receive 
contributions from an infinite series of multi-loop diagrams, with a specific structure 
(essentially, chains of ladder diagrams). Recently, the first few such diagrams for the 
polarization tensor have been explicitly computed by Bodeker ^2|, with results which 



agree with the (first order iteration of the) solution to the Boltzmann equation. But it 



p'=p+q 



k k'=k-q 

Figure 1: Elastic scattering in the (resummed) Born approximation. The continuous hues 
refer to hard gluons (these are off-equihbrium propagators), while the wavy line is the 
soft gluon exchanged in the collision. The blob stands for HTL resummation. 

is clear that, in general, computing directly these diagrams would be a tedious exercice, 
especially since important cancellations occur among the various graphs |^ ^; this will be 
further discussed in Sees. 3.2 and 4.1 below, in) The resummation of the collision effects 
drastically modify the longwavelength behaviour of the transverse colour conductivity 
ari^^ = 0,p — > 0): the would-be divergence of the HTL result for ax, namely a^ oc mj^/p, 
is now screened away by 7, with the net result that axiuj = 0,p — ^ 0) = m|,/3(7 — 5). Here 
(5 is a term of order g'^T (and which satisfies \6\ < 7), to be computed in Sec. 4.2 below. 
iv) For ultrasoft momenta p ~ g'^T, the above formula for ax holds up to corrections of 
0(ln"^) (with In = ln(l/c/)). 

2 The Boltzmann equation 

In this section, we review the main features of the Boltzmann equation for colour relax- 
ation which was recently derived in Ref. 0]. There are no new results to be reported 
here, but the equations derived in 0] will be presented in a slightly different way, to better 
emphasize the difference between coloured and colourless excitations, and also between 
Abelian and non-Abelian gauge theories. Moreover, the diagrammatic representation of 
the collision term will be explained in more detail, in order to facilitate the discussion of 
the diagrammatic interpretation of the ultrasoft amplitudes, in Sec. 3.2 below. 

The space-time inhomogeneities in the distribution of the hard particles (transverse 
gluons) are described by a density matrix 6Gab{k,x) where the momentum k is hard 
{k ~ T) and on-shell (fcg = k^), and the derivative dx is ultrasoft {dx ~ g^T). Below, we 
shall be interested either in colourless fluctuations (in which case 6Gab = ^ab^G), or in 



coloured ones in the adjoint representation (such as 5Gab = if^'^SGc)- Also, we shall find 
convenient to use the following parametrization for the density matrix, where the on-shell 
structure is explicit: 

6Gab{k,x) = -po{k)Wab{k,x)—— 

dfco 

= Ppoik)Wab{k,x)No{ko)[l + No{ko)]. (2.1) 

In this equation, A^o(^o) = l/(e^'^° — 1) and po{k) = 2TTe{ko)6{k'^) are, respectively, the 
thermal distribution and the spectral density for hard transverse gluons, and the new 
function Wab{k,x) has support only at the mass-shell: Wab{k,x) = 9{ko)Wab()<;X) + 
9{—ko)Wba{—y^, x). [The density matrix 6Nab(k, x) in the Introduction, eq. (|1.3|) , is related 
to the function Wab{Kx) by 6Nab()<.,x) = -WabO<;X){dNo/dk).] 

The Boltzmann equation is the kinetic equation satisfied by the density matrix to 
leading order in g [Q. It reads (in matrix notations): 

k-D^, 6G{k,x)] - 2gk^'F^^{x)dlG^{k) = G{k,x). (2.2) 

In the l.h.s., k ■ Dx is the gauge- covariant drift operator, with D^ = 5^ + igA^ and 
(9^ ~ ^A ~ g'^T, so that D^ = 0{g^T); k^'F^^{x)^^, with F^^ = [D^, D^]/ig, is the "force" 
term acting on the equilibrium correlation function (^^(A;): 

G<ik) ^ poik)Noiko), G>{k) ^ p,{k)[l + iVo(fco)]. (2.3) 

(The second function G^lk) will be needed below.) 

In the r.h.s. of eq. ( |2.2| ), G{k,x) is the collision term associated to the one-gluon 
exchange process depicted in Fig. 0. A priori, all the lines in this figure (that is, both 
the external lines associated with the colliding particles, and the wavy line associated to 
the exchanged gluon) are off-equilibrium propagators. However, to the order of interest, 
the collision term can be linearized with respect to the off-equilibrium fluctuations in the 
propagators of the external lines, and the internal propagator can be taken to be the 
equilibrium propagator. Since the collision term for colour relaxation is dominated by 
soft momentum transfers {g^T ^ q ^ gT) ||, ||, ||, |^, ^, the propagator of the exchanged 
gluon has to be dressed with the corresponding hard thermal loop |16|, |17 . 



The scattering process in Fig. |^ can be associated to the coUisional self-energy in 
Fig. 1^ (see Ref. 0] for more details on the formalism). Upon linearization, this leads to 
the four processes displayed in Fig. ^. Each diagram involves fluctuations in one of the 
four external lines in Fig. [l|. Thus, G = Ci + G2 + C^ + C4, where Giik^x) involves the 
fluctuations 6G{k,x) in the incoming fleld with momentum k (Fig. ^a), and G2, C3 and 
G4 involve fluctuations along the lines with momenta k', p and p' (Figs. |^.b, c and d, 




k k-q 

Figure 2: Self-energy describing collisions in the (resummed) Born approximation. All 
the lines represent off-equilibrium propagators. The continuous lines refer to the hard 
colliding particles in Fig. 0. The wavy lines with a blob denote soft gluon propagators 
dressed by the screening effects. 

respectively). In these figures, the off-equilibrium propagators are marked with a cross; all 
the other lines denote equilibrium propagators. In particular, Ci{k,x) = —T{k) SG{k,x), 
where 



r(A;) = S< (fc) - S> (fc) = -2ImS^(A;) 



eq\ 



(2.4) 



is the quantity which determines the quasiparticle damping rate 7 = r(/i;o = k)/{Ak) 
g^T\n{l/g) p|,|0|,|T],|2l. 



By using the parametrization ( p.lD for the density matrix, the (linearized) collision 
term can be compactly written a^: 



(-^abyk, X) 



dT \M^k^M^ No{ko) Noipo) [1 + No{k',)] [1 + iVo(p[,)] 
x{N[NWab{k,x) - {T^T'),aWcd{k',x)) + 
+ {T^T\,{rr),s{Ws,{p,x) - W,s{p',x))}. 



(2.5) 



In this equation, |A^pfc^p'fc'| is the matrix element squared corresponding to the one- 
gluon exchange depicted in Fig. [^, and dT is a compact notation for the measure of the 
phase-space integral: 



dT = f3 



d'^p f d'^q 



{2ny J (27r) 



Poik)poip)poip + q)poik - q). 



(2.6) 



The four terms within the braces in eq. ( p. 51) are in one to one correspondance with the 
diagrams |3.a, b, c and d. The appearance of the matrix element squared, and also of 



•''Eq. (2.5) is merely a convenient rewriting of eqs. (3.99)-(3.104) in Ref. 



■^ 




(a) 



k-q 




(b) 



>^ 




k-q 



k-q 



(C) 



(d) 



Figure 3: Pictorial representation of the linearized collision term. Each one of the four 
diagrams correspond to off-equilibrium fluctuations in one of the colliding fields (the one 
which is marked with a cross). All the unmarked propagators are in equilibrium. 



the various equilibrium statistical factors in eq. ( |2.5|) , is familiar. What is specific to the 
problem at hand is the colour structure in eq. ( p.5|) , which is at the origin of an important 
difference between coloured and colourless excitations: 



Consider first the case of a colourless fluctuation, for which 5G 



ab 



Wab = SabW. The various colour traces in eq. ( WM are trivial (e.g., (T"-T^ 



SabSG, and 
= N6ab), so 



that C, 



ab 



6nhC, with 



C{k,x) = -N^ J dr\Mpk^p>k'f No{ko)No{po)[l + No{k'o)][l + No{p',)] 
X [W{k, x) - W{k', x) + W{p, x) - W{p, x)}. 



(2.7) 



This is the standard collision term for one-gluon exchange used in previous applications 
of kinetic theory to the hot quark-gluon plasma or to the electroweak plasma [ |12| |. 

What is remarkable about eq. ( pIT] ) is that the corresponding phase-space integral 
is dominated by relatively hard momentum transfers gT ^ q ^ T, even though each of 
the four individual terms in the r.h.s. is actually saturated by soft momenta. This is a 
consequence of the cancellation of the leading infrared contributions among the various 
terms [Q. For instance, for soft q, W{k',x) = W{k — q,x) ~ W{k,x), so that the IR 
contributions to the first two terms in eq. ( |2.7| ) cancel each other. This corresponds to 
a cancellation among the graphs displayed in Figs, ^a and b, to be further discussed in 
Sec. 3.2 below. A similar cancellation occurs between the last two terms in eq. 



8 



namely W(j),x) and W{p',x). Thus, in order to see the leading IR (g ^ T) behaviour of 
the full integrand in eq. ( p.5|) , one has to expand W{k',x) and W{p',x) to higher orders 
in q. This generates extra factors of q which remove the most severe IR divergences in the 
collision integral. (This is the familiar (1 — cos^) factor of the "transport cross section" 



33[|.) As a result, the typical rate involved in the calculation of the transport coefficients 
like the shear viscosity is ~ g^Tln{l/g), where the logarithm originates from screening 
effects at the scale gT |^, |ll|. This is suppressed by one power of a = g^/4:7i with respect 
to the damping rate 7 ~ g'^T\n{l/g). 

Consider now the case of colour fluctuations corresponding to a density matrix 
W{k,x) of the form W{k,x) = Wa{k,x)T°-. The colour algebra in eq. (|2.5|) can be 
performed with the following identities: 

TV _ /V^ 

Tr(T'^T^T^) = zr'^ - , {T'^T'UT^T^UiTndd = ^/"'^ -j- ■ (2.8) 

The resulting collision term is of the form C = CaT"" with 

Ca{k,x) = -N^ j dT\M,k^,yk'? N,{k,)N,{p,)[l + N,{k',)][l + N^{p'^)] 

yi^Wa{k,x)- ^Wa{k',x)- ^{Wa{p,x) + Wa{p',x))Y (2.9) 



There are two notable differences with respect to eq. (|2.71 ): 

i) The flrst two terms within the braces enter with a relative factor 1/2, so they do not 

cancel each other when g -^ : 

Wa{k,x)- ^Waik-q,x) - ^WaiKx). (2.10) 

Rather, their overall contribution is half the corresponding contribution of the flrst term 
alone, that is, r{k)/2. 

ii) The last two terms in eq. ( p.9|) enter with a factor 1/4 and add each other. This is so 
because the colour matrix W^^{p, x) in eq. ( p.5|) is antisymmetric, rather than symmetric, 
as it is for colourless fluctuations. Accordingly, for soft q, 

^(Wa{p,x) + Wa{p + q,x)) ^ ^Wa{p,x). (2.11) 

Thus, for colour fluctuations, the colour structure of the collision term prevents a complete 
cancellation of the leading infrared contributions: like the damping rate, the collision term 
for colour relaxation is saturated by soft momentum transfers {g^T ^ q ^ gT), for which 
eqs. ( p.lO|) and (|2.11| ) hold and the collision term (pl9| ) simplifles to 0, 0: 



N^ r ,^. . , ,9 dA^n dA", 



a(fc,x)^-— f dT\M,k^,,,,f —^^ {Wa{k,x) - W,{p,x)}. (2.12) 

2 J akn d»n 



In the same approximation, the matrix element |7Wp can be evaluated as: 

\M\' = IQg'elel \*Vi{q) + (q x v) • (q x v') *V,{q)\\ (2.13) 

where v = k, v' = p, and *T>i and *Vt are the longitudinal (or electric) and the transverse 
(or magnetic) components of the (retarded) gluon propagator, in the hard thermal loop 
approximation [|l^, |l^]. The phase-space measure ( p.6|) can be similarly simplified. This 
eventually yields a simpler equation for the density matrix VFa(k, x) which, remarkably, 
is consistent with VFa(k, x) being independent of the magnitude k = |k| of the hard 
momentum. That is, 

Waik,x) = gWaix,v), (2.14) 

where v = k/fc is the velocity of the termal particle (a unit vector), and a factor of g has 
been introduced to keep in line with the notations of Ref . [Q . 

Finally, the Boltzmann equation, written as an equation for Wa(a;, v), reads [Q: 

{v D,r'W,{x,v) = vE'^(x)-m^^^y'^$(vv'){iy"(x,v)-W^"(x,v')}. 

(2.15) 

The angular integral above runs over all the directions of the unit vector v', and mj^ is 
the Debye mass squared: 

rrijj = — / dpp — — = — - — . (2.16) 



TT^ Jo dp 

Furthermore: 

$(vv') = (27r)Y^%o-q-v)%o-q-vO|*2^z(g) + (v,-vO*A(g)f, (2.17) 

with the two delta functions expressing the energy conservation at the two vertices of 
the scattering process in Fig. |l]. Up to a normalization, the function ^(v ■ v') represents 
the cross section for the collision between two hard particles with velocities v and v' 
exchanging (in the t-channel) a soft (dressed) gluon. 

The collision term in eq. ( |2.15| ) involves two pieces: one which is local in v (pro- 
portional to Vr"(x, v)), and one which is non-local (involving the kernel $(v • v'). The 
coefficient of the local piece is proportional to F : 

2 g^NT fdn', . T(ko = k) 
ml) / $ V ■ V = ^ = 7. 2.18 

By using the expression above, we can rewrite the Boltzmann equation (|2.15|) in the 
following way: 

{v D^Y^Wb{x,v) = vE"(x)-7{iy"(x,v) - (iy"(x,v))}, (2.19) 

10 



which emphasizes the fact that the quasiparticle damping rate 7 sets the time scale for 
colour relaxation: t^oi ~ I/7 ~ 1 / {g^T \ti{1 / g)) 0] (see also Sec. 4.2 below). In eq. (|2.19| ) 



we have introduced a notation which will be used hereafter: for an arbitrary function of 
V, say -F(v), we denote by (-F(v)) its angular average with weight function $(v • v'): 

/ f^1>(vvOF(v O 
which is still a function of v. 



in^)) ^ -^finr^JZT^TJT^ ^ (2-20) 



We conclude this section by recalling that eq. ( |2.15| ) is invariant under the gauge 



transformations of the background field, and also with respect to the choice of a gauge 
for the shortwavelength fluctuations (here, the hard {k ~ T) fields which take part in the 
collective motion and the soft {g^T ^ q ^ gT) gluons which are exchanged in the collision 
process). In Ref. [Q, eq. ( p.l5| ) was derived in Coulomb gauge, but we expect it to be 
gauge-fixing independent. Except for the collision term, this has been explicitly verified 
in |jl|] (see also Refs. |T^, 0). The collision term should be gauge-fixing independent 
as well, since it involves only the off-equilibrium fluctuations of the (hard) transverse 
gluons, together with the (gauge- independent) matrix element squared ( ^.13| ). However, 
an explicit proof comparable to the corresponding one for the non-Abelian Vlasov equation 
|l[] is somewhat tedious: in an arbitrary gauge (e.g., a covariant one), one has to consider 
collisions involving fictitious degrees of freedom (hard longitudinal gluons and ghosts), 
and verify that their respective contributions to the collision term mutually cancell. 



3 Ultrasoft amplitudes 

In this section, we introduce and study the ultrasoft amplitudes, i.e., the contributions 
to the one-particle irreducible amplitudes with ultrasoft external lines which are obtained 
from the solution to the Boltzmann equation. 

3.1 The induced current 

The longwavelength colour fluctuations of the hard particles generate a colour current 
given by (the factor of 2 below stands for the two transverse polarizations): 

r,{x)=2gJ-0^^k,TT(T^6G{k,x)), (3.1) 

which acts as a source in the Yang-Mills equations for the ultrasoft colour fields A^: 

(D^F^.Tix) = j;{x). (3.2) 
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By using the parametrization ( |2.1| ) for the density matrix, one can perform the integral 



over the radial momentum k = |k| to obtain: 

j^iix) = mlJ^v'^W^{x,^r), (3.3) 

with the Debye mass rriu defined in eq. ( ^^.161 ). By using the equation of motion ( |2.15| ) 
for iy"(x, v), one can verify that the current ( |3.3| ) is covariantly conserved, 

D,f = 0, (3.4) 

as necessary for the consistency of the mean field equations of motion (|3.2|) (recall that 
Di'D'^F^^ = 0). Indeed, eq. {^J§ implies: 



D,f = mlj — wW{x 



g'^NT fdn fdn' 



^ 2 



/ £;^ / ^ $(v ■ v'){iy"(x, v) - W\x, v')}, (3.5) 

which is zero because both terms in the r.h.s. vanish after the angular integration. 

By solving the Boltzmann equation, one can obtain the density matrix Wa{x,\), 
and therefore also the induced current Ja{x), as functional of the gauge fields A'^{x). 
Since eq. (p.l5|) is non-linear with respect to the fields A'^, the resulting functional Ja[A] 



j- = njA;; + -rf;^A,M^+... (3.6) 



will be non-linear as well, and can be formally expanded as follows: 

1 

2 

The coefficients in this expansion are the one-particle-irreducible amplitudes of the fields 
A^^, evaluated in thermal equilibrium |l^, |^, |l|]. For instance, 11^^ = S^-^U^i, is the polariza- 
tion tensor, F^^*^^ = if°''^'^T ^i,p is a correction to the 3-gluon vertex, etc. These amplitudes 
will be referred to as the ultrasoft amplitudes. 

Some of the properties of the ultrasoft amplitudes follow immediately from the 
previous discussion: For generic momenta of order g^T., they are of the same order in 
g as the hard thermal loops [|^, |19|, |l|, |1^, |1^ |[, which they generalize by including 



the effects of the collisions. Furthermore, they are gauge-fixing independent (like the 
Boltzmann equation itself), indicating that only the physical hard degrees of freedom 
of the plasma (namely, the on-shell transverse gluons) contribute to these amplitudes. 
Also, they satisfy simple Ward identities which follow from the conservation law (|3.4|) by 
successive differentiations with respect to the fields A'^. For instance: 

p^n^,(p) = 0, 

P^T^,p{P,Q,R) = n,p(Q) - n,p(i?) . (3.7) 

All these properties are, of course, very reminiscent of the hard thermal loops, and, as we 
shall see later, there are other similarities. But let us first discuss the interpretation of 
the ultrasoft amplitudes in terms of Feynman diagrams. 
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3.2 Diagrammatic interpretation of the ultrasoft amplitudes 

In this subsection, we discuss the interpretation of the solution to the Boltzmann equation 



in terms of Feynman diagrams. (See also Refs. |3^, ^ for a related analysis in the context 
of scalar field theory, and Ref. |3^ for a recent calculation of some of the diagrams relevant 
to QCD, namely those in Fig. |^ below.) This analysis will show that, unlike the HTL's 
— which correspond to one-loop diagrams |l^, |19| — , the ultrasoft amplitudes receive 
contributions from an infinite set of multi-loop Feynman graphs, which, in the kinematical 
regime of interest, contribute all at the same order in g. 

Our discussion here will be only qualitative: we shall not compute Feynman graphs 
explicitly, but rather rely on the diagrammatic representation of the collision term (cf. 
Figs. 0, and ^ in order to identify, by iterations, the structure of the diagrams con- 
tributing to the ultrasoft amplitudes. 

To carry out the analysis, it is convenient to use the original form of the Boltzmann 



equation, where the collision term is directly related to the self-energy. This is eq. (2.2) 
with the collision term (|2.5| ), which we rewrite here as follows: 



vD,, 6G{k,x)\ = gv^F^,{x)d'^G<{k) + C{k,x), 

C = Ci + C2 + Cs + C4, Ci = - -^ r(A;) 6G{k, x), (3.8) 

where v^^ = k^ /k^, and the four pieces of the collision term correspond, respectively, 
to the linearized fluctuations depicted in Figs. |^.a, b, c and d. (Note that the present 
normalization of the collision term differs by a factor l/2A;o from the previous one in 
eq. (|2.2|) .) For comparison with perturbation theory, it is useful to regard the collision term 
as a "small perturbation" and solve the Boltzmann equation (|3.8|) formally by iterations. 



The zeroth order iteration is the solution to eq. ( |3.8| ) with the collision terms ex- 
cluded: 

5G^''^ = g^v'^F^,d''G<. (3.9) 

V ■ U 

Here, l/(f ■ -D) is a compact, but formal, notation for the retarded Green's function 

A/j(x, y; v) of the covariant drift operator v ■ D. This satisfies: 

{v D^U/\'^{x,y-Y) = <5-<5(4)(x-i/), (3.10) 

with A/j(x, y; v) = for a;o < yo, and has the following expression (with t = xq — yo): 

A?*(x,i/;v) = e{t)5^'H^-y-vt)U'^\x,y) ^ Ua—^-—-y,b), (3.11) 

where f/(x, y) is the Wilson line connecting the points x and y : 

U{x,y) = e-''S'^''^-^'\ (3.12) 
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Figure 4: Mean field approximation, or zeroth order iteration for the solution to the 
Boltzmann equation, eq. (pl9|). 



and the integration path in eq. ( |3.12|) is fixed by the delta function in eq. ( p.ll| ). Aji{x, y; v) 
is the eikonal propagator along the straightline trajectory of velocity v. 

Eq. ( p.9| ) can be given the diagrammatic representation in Fig. ^ where, for more 
clarity, we have distinguished the insertion of the electric mean field E^, the "Lorentz 
force" in the r.h.s. of eq. (|3.8| ), from the insertions of colour fields A'^ due to the covariant 
derivative v ■ D in the l.h.s. Thus, the propagator on the left of the electric field E^ 
is the eikonal propagator (|3.11|) , while the propagator on the right is dkGQ^k). The 
corresponding colour current, namely: 



^^''{x) 



29 



/7Si^'^Tr(T«5G(0)(^,x))=rn?,/^ 



dn v^v'E' 



V D 



(3.13) 



involves a supplementary integration over the hard momenta /c, which, in terms of di- 
agrams, corresponds to closing the straight line in Fig. ^ into a hard loop. Thus, the 
polarisation amplitudes generated by j^°^ (cf. eq. (|3.(j|)) are one-loop amplitudes where 
the internal momentum is hard, while all the external lines are soft (or ultrasoft); some 
examples are shown in Fig. ^. These are precisely the hard thermal loops, which have been 

n)|. As well known, the HTL's 



originally computed from one-loop diagrams indeed 
are only a part of the corresponding one-loop amplitudes |18|, |l^, [l^ (namely, the leading 
order part for soft external lines), and this part is directly singled out by the coUisionless 
kinetic equation, cf . eqs. ( |3l9| ) and ( |3.13| ) |1|, [1^, ^ . 



Consider now the first order iteration of the collision term, which yields: 



5G(i) 



1 



V D 



C[(5G(°)] . 



(3.14) 



The collision term C[(5(j'-'^^] is illustrated in Fig. ^, which should be compared to Fig. ^. 
For simplicity, we have only represented here single field insertions, that is, we have lin- 
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Figure 5: Two and four-gluon vertices in the HTL approximation. 
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Figure 6: The first order iteration of the colhsion term C[5G''''^^], to hnear order in the 
background field. 
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(c) (d) 

Figure 7: First order iteration for the ultrasoft polarization tensor. 

earized C[(5(j^°^] with respect to the colour mean field. This is all what we need in order 
to compute the first-order iteration of the ultrasoft polarization tensor H^u. The corre- 
sponding result is illustrated in Fig. ^ and involves loop corrections to the corresponding 
HTL (cf. Fig. |.a). 

Before going on with higher iterations, let us make some comments on the diagrams 
in Fig. 0. These should be regarded as diagrams of the thermal perturbation theory in real 
time, or linear combinations of them. For instance, the self-energy insertion in Fig. ^a 
stands for the combination r(A;) = T,fg{k) — T,^q{k) (cf. eq. ( |2.4| )) and thus corresponds to 
the insertion of the quasiparticle damping rate in the hard internal line. Similarly, the soft 
internal fine (with a bubble) in Fig. |^.b stands for either *'D^^{q) = -O^VfJl^f^^*!) 

or *V>Xq) = 



flU 



HV^ 



\0) 



*P/jnfo.*pA , where Wf^. denotes the two-point HTL (cf. eq. (|]T 



fJ,W 



below), and *V^iy is the HTL-resummed propagator (the subscripts R and A refer, as usual, 
to retarded and advanced propagators) [||. To simplify the graphical representation, it 
is convenient to replace these graphs with the corresponding ones in the imaginary time 
formalism, where the two diagrams in Figs. |^.a and b are replaced by the graphs in Fig. ^a 
and b, respectively, while the diagrams in Figs, ^c and d remain formally the same, and 
are globally represented in Fig. |^.c. 

Now, since the Boltzmann equation has been obtained from the exact field equations 
by using various kinematical approximations [Q, the correspondence between its solution 
(here, eq. (|3.14|) ) and the diagrams in Fig. ^ is only approximate: the kinetic equation 
isolates only the dominant parts of these diagrams for ultrasoft external lines. In fact, in 
the same way as the Vlasov equation ( p.9|) provides the leading contribution to the one 
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(b) 




(c) 



Figure 8: A simpler redrawing of the diagrams in Fig. ^ 



loop diagrams when the external momenta are ~ gT, the Boltzmann equation generates 
automatically the leading contributions to the ultrasoft amplitudes when the external 
momenta are ~ g^T. Recently, this has been verified explicitly by Bodeker |^2[, who 
computed the diagrams in Fig. ^ and got the same result (namely, eq. ([4.19|) below) as 
that obtained from the iteration of the Boltzmann equation. In fact, the approximations 
leading to the Boltzmann equation |Q and those performed in Ref. |^ are similar. They 
all rely upon the following chain of inequalities: 



c^j; < g < fc. 



(3.15) 



which are controlled either by powers, or, at least, by a logarithm of the coupling constant. 
Specifically: 

(a) The gauge covariant gradient expansion retains the terms of leading order in 
dx/k] this is an excellent approximation since the neglected terms are of 0((7^) or less. 
Diagramatically, this translates into the fact that the smooth lines in Figs. ^ |^, or ^ (and 
also in the diagrams to come) represent ezA;ona/ propagators, of the form (cf. eq. (|3.11| )): 



A«(P,v) 



(3.16) 



V ■ P + ie 

rather than standard tree-level propagators. Moreover, all the vertices in these diagrams 
are simplified by systematically ignoring the external momentum P. 

(6) In the construction of the collision term, we have retained only the terms of 
leading order in an expansion in powers of dx/q. This approximation, which assumes that 
both particles taking part in the collision (see Fig. |1|) feel the same mean field, is needed 
in order to put the collision term into a form local in x. But for colour fiuctuations at 
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the scale g'^T, this is correct only up to corrections of 0(l/ln(l/5f)) : indeed, dx ~ g'^T, 
while the cross section in eq. ( p.6|) is logarithmically sensitive to momenta q ~ g^T (see 
eq. ( |3.29| ) below). Diagrammatically, this affects only the diagram ^.c (more generally, 
the diagrams involving two or more hard loops; see, e.g.. Fig. p!0| below) , where it amounts 
to assume that both the internal wavy lines (the two gluons connecting two hard bubbles) 
carry the same momentum q. (Strictly speaking, if one of these lines has a momentum q, 
then the other one should rather carry a momentum p — q.) 

(c) Within the collision term, we have neglected, wherever possible, the exchanged 
momentum q as compared to the hard momenta of the colliding particles (recall the 
discussion after eq. ( p.9| )). This is a good approximation since q/k = 0{g) or less. Dia- 
grammatically, this entails more simplifications in the propagators and vertices in Figs. |: 
the velocity remains unchanged when running along a given hard loop (e.g., in Fig. |^.c, 
there are only two velocities: v = k for the left hand loop, and v' = p for the right hand 
one), and the momentum q is neglected in all the vertices. 

It is interesting to examine the validity of these approximations in the separate 
cases of colour fluctuations and colourless ones. The approximations (a) and (6) are 
quite generic in relation with the Boltzmann equation, and are actually better justified 
for colourless fluctuations than for coloured ones: Indeed, we have seen in Sec. 1.1 
that the colourless fluctuations relax mainly via hard (or large angle) scattering, g ~ T, 
with a typical relaxation rate ~ g'^Tlnll/g). In this case, the effect of the collisions 
becomes a leading order effect only for very soft inhomogeneities, dx ~ g^T, for which 
both inequalities dx <^ k and dx '^ q are very well satisfied. On the other hand, the third 
approximation (c) does not apply to colourless fluctuations, for which, because of the 
cancellations discussed after eq. (|2^ ), g ~ fc ~ T. Note that it is precisely approximation 
(c) which allowed us to reduce the original collision term ( p.5| ) to the simpler expression 
in eq. ( ^121) . 

The simplifications arising from approximation (c) (cf. eqs. ( |2.10D and ( p.ll|) ) have 
actually a simple diagrammatic interpretation as cancellations among Feynman graphs. 
The two diagrams in Figs. |^.a and b correspond respectively to the terms involving W{k, x) 
and W{k',x) in eq. ( p^.5| ). For colourless fluctuations in QCD, or, equivalently, for electric 



fluctuations in QED, eq. ( p.7|) teaches us that these two diagrams cancel each other in the 
limit where q is neglected next to k or p. That is, each of these diagrams is individually 
dominated by soft momenta q, but their leading infrared contributions mutually cancel 
in the sum of the diagrams, so that we are left with the (subleading) contribution of hard 
q momenta. (This is what makes the colourless collision term ( p.7|) difficult to deal with; 
see, e.g., p|,|,[ni.) 

In QED, this cancellation has been also verified via direct diagrammatic calculations, 
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(a) 



(b) 



Figure 9: (a) A ladder diagram, as generated by iterations of the first two pieces, Ci and 
C2, of the colhsion term; the smooth hnes with a bubble are eikonal propagators dressed 
with a damping rate 27. (b) The sum of all the ladders in (a), as obtained after using the 
partial cancellation between vertex and self-energy corrections to effectively remove the 
ladders; the thick line is an eikonal propagator with a damping rate 7. 



in Refs. ||2^, ^, ^, ^]. In particular, in Ref. |^, this has been related to the absence 
of HTL vertices with four external photons. Indeed, the two diagrams in Figs. |^.a and 
b can be generated from the four-particle hard thermal loop in Fig. |^.b, by closing two 
of the external lines in all the possible ways. In QED, this involves the four-photon 
HTL which, however, is well known to vanish |l^, |ll : the HTL-like contributions of 
the individual diagrams with four external photons (as in Fig. ^.b) mutually cancel after 
summing over the permutations of the external lines. In the present framework, the sum 
over the permutations corresponds precisely to the sum of the two diagrams in Figs. |[a 
and b, so this sum has to vanish as well. 

In QCD, on the other hand, the sum over permutations produces a colour commu- 
tator, which thus provides both a non- vanishing four-gluon HTL |]TB|, |TU|, and a non-zero 
global contribution from the diagrams in Figs. |^. This is the content of eqs. ( |2.9| )-( |2rT2| ). 
In fact, eq. (|2.10| ) shows that, even in QCD, there remains a partial compensation be- 
tween the self-energy and vertex corrections in Figs, [^.a and b (while the two diagrams in 
Fig. ^c and d rather reinforce each other; cf. eq. ( p. 11 )). The net effect is that /ia// of the 
self-energy correction is cancelled by the vertex correction, with the factor 1/2 coming 
from the colour algebra (cf. the first trace identity in eq. (|2.8|)). 



We now turn to the higher order iterations. Clearly, by iterating the self-energy 
insertion in Fig. |[a, one ends up with replacing the propagator of the hard gluon with a 
dressed propagator which includes the damping rate. That is, the bare eikonal propagator 
( |3.16D is replaced by the following dressed propagator 



*A^,(P,v) 



t; ■ P + 2i7 ' 
to be graphically represented by a straight line with a blob (see Figs. 



(3.17) 



10). 



|.a and 

Equivalently, this resummation can be achieved by moving the first collision term Ci into 
the l.h.s. of the Boltzmann equation (p.8|). Similarly, by iterating the vertex correction in 
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Figure 10: A generic ladder diagram contributing to the ultrasoft polarization tensor, as 
obtained from the Boltzmann equation. 

Fig. |^.b one generates the ladders diagrams depicted in Fig. ^.a. Finally, diagrams with 
two or more hard loops will be generated by iterating the other two pieces, C3 and C4, of 
the collision term (cf. Figs. 0.c and d, and Fig. ^.c). 

We conclude that the typical diagrams which are resummed by the solution of the 
Boltzmann equation ( p.8|) are as shown in Fig. |10|. They involve a chain of an arbitrary 
number of hard loops, each of them dressed by ladders and damping effects as in Fig. |^.a, 
and connected one to the other by pairs of soft gluons. The smooth lines with a blob 
represent the dressed eikonal propagator (|3 . 1 7|) , while those without a blob are thermal 
correlation functions like G^ and Gq (cf. eq. ( |2.3| )), or derivatives of them (cf. eq. ( |3.9| )). 
Now, all the previous examples involve diagrams which contribute to the ultrasoft polar- 
ization tensor (or 2-point amplitude) Il^^{P). But, of course, similar diagrams exist for all 
the higher point ultrasoft vertices: they can be obtained by inserting more external lines 



along the hard loops in Fig. |T0|, on any of the internal eikonal lines (which, in general, 
can be seen as eikonal propagators in a background field; cf. eq. ( |3.11|) ). 

It is finally possible to give a simple graphical interpretation of the (partial) infrared 
cancellations between self-energy and vertex corrections, as discussed above. To the order 
of interest, the only effect of the ladder corrections in Fig. ^.a is to reduce to damping 
rate0 in eikonal propagators like ( ^.171 ) from 27 to 7 (cf. eq. (^TTD|)). This is depicted in 
Fig. ^.b, where the thick internal line denotes the following eikonal propagator (compare 
to eqs. (|3l6|) and (^TzD): 



Ai?(P,v) 
while the thin line corresponds to 9fcG^(A;). 



V ■ P + ij 



(3.18) 



3.3 General properties and iterative solutions 

We now return to a discussion of the general properties of the ultrasoft amplitudes. Since 
they encompass, and generalize, the HTL's, we expect these amplitudes to describe phe- 



^In QED, the sum of all the self-energy and ladder corrections depicted in Fig. |^.a simply vanishes to 
the order of interest Q ; cf. the discussion after eq. (^.7D . 
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noinena like Debye screening or Landau damping (possibly modified by the effects of the 
collisions), and also transport phenomena, which are made possible by the collision term; 
the exemple of the colour conductivity will be discussed in the next section. 

In order to look for Debye screening, it is enough to consider static fields, that 
is, colour field configurations which are described by time-independent vector potentials 
A^(x). In this case, the ultrasoft amplitudes reduce to the usual Debye mass term m|, = 
g'^NT'^/3 for electrostatic fields, as obtained from the HTL's. In order to see this, it is 
convenient to decompose the functions W°'{x,\) as follows [Q : 

W%x,^r) = -A^oi^) + A''{x,^r). (3.19) 



From qs. ( |2.19| ) and ( |3.19| ), the following equation is obtained (recall that E^^ = -Da6^6 ~ 
d'Al) : 

{vD,r'Ab{x,v) = d'ivA^) - i{A''{x,v) - (^'^(x,v))}, (3.20) 

while eq. ( |3.3| ) shows that the current can be rewritten as: 

jiiix) = -6^'mlA^,{x) + mlJ^v'^Aa{x,^r). (3.21) 

In obtaining eq. (p.20 ), we have also used the fact that: 

Al{x)-{Al{x)) = 0, (3.22) 

since the collision term vanishes for any function which is independent of v. 

In eq. ( p.20|) , the time derivative of the vector potentials (i.e., the term d^{v-A)) acts 



as a source for the functions A°'{x, v). Since we are looking here for solutions which vanish 
in the absence of sources, it follows that A°'{x, v) = (and therefore W°'{'x, v) = — ylg(x)) 
when the gauge potentials are time- independent. Then, eq. ( |3.21|) reduces to: 



j^(x) = -V"^D^o(x), (3.23) 

which is the same expression as in the HTL approximation |T| . That is, for static external 
legs, all the ultrasoft vertices with n > 3 external lines vanish, while Ilfj,u{uj = 0,p) = 

Eq. ( p.23| ) shows, in particular, that the value of the Debye mass is not modified by 
the collisions among the hard particles. An alternative derivation of this result has been 
recently given in Ref. ^^ (see also Sec. 4.2 below). This is not unexpected since we know 
^ that the first correction to mj^, of 0{g^T'^ln{l/g)), comes out from the interactions 



between soft and ultrasoft fields. 
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For time-dependent fields, however, the colhsions among the hard particles do play 
a role, and, at very soft momenta P ^ g'^T (by which we mean that both the frequency 
P*^ = u;, and the spatial momentum p = |p|, are of order g'^T or less), they can even 
dominate over the mean field effects. This may be seen by considering the formal solution 
of the Boltzmann equation (|2.19| ) obtained by iterations. There are several ways to 



organize the iteration. For instance, we may iterate the whole collison term in the r.h.s. of 
eq. ( |2.19|) , similarly to what we have done in the previous subsection. Since the collision 



term is proportional to 7, the resulting solution is a formal expansion in powers of 7. 
Specifically, we write W = W^^^ + W^^'^ + W^^'^ + . . ., where W^°'^{x,v) satisfies the 
transport equation in the mean field approximation (or Vlasov equation) 



\abTjAO) 



{vD^rWl''>{x,v) = vE^{x), (3.24) 

while the A^th order correction Wjj^^\x, v) is proportional to 7^. The (retarded) solution 
to eq. (|3.24| ) involves the eikonal propagator Aij(x,y;v), as defined in eq. ( p.ll| ). Thus, 




in compact notations: 

W^(i) = 

This expansion maintains explicit gauge symmetry at each order in 7: indeed, since 
both pieces of the collision term (i.e. the local piece —'jW°-{x,\) and the non-local one 
'y{W°'{x, v))) are treated on the same footing, the current conservation law ( |3.4| ) is verified 
at each step in this iteration. Then, e.g., the polarization tensor Iljf^'' constructed in the 
A^th iteration is guaranteed to be transverse. 

The discussion in the previous subsection provides us with the diagramatic inter- 
pretation of the expansion (|3.25 ). The zeroth order solution ty(°) corresponds obviously 



to the mean field, or HTL, approximation (cf. eqs. ( |3.9| ) and ( |3.13| ), and Figs. ^ and ^. 
The first order iteration W^^^ corresponds to the diagrams in Figs. or |[ Specifically, 
the expression of W^^^ in eq. ( p.25| ) involves two pieces within the braces in its r.h.s.: the 
first piece corresponds to the sum of the self-energy and vertex corrections depicted in 
Figs. |^.a and b (or, equivalently, in Figs. |.a and b); similarly, the second piece in W^^^ 
corresponds to the sum of the two diagrams in Figs. 0.c and d. But it is only the set of the 
four diagrams in Fig. |^ which is globally gauge invariant and provides a transverse contri- 
bution U^^} to the polarization tensor [^ (see also Sec. 4.1 below, especially eq. ( [4.19] )). 
Note that, for ultrasoft gradients D^ ~ g^T ~ 7, the expansion ( p.25| ) is a formal one: 
indeed, all the terms are equally important. 
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A different expansion is obtained by choosing only the second piece of the colhsion 
term, namely 'y{W°'{x, v)), as the perturbation. To do this, we move the term —'yW°'{x, v) 
to the l.h.s. of the Boltzmann equation, and define the "dressed" eikonal propagator: 

Af{x, y; v) = e{t) d'^'^ (x - y - vt) e"^* U^\x, y) ^ {x, a ' . y, b). (3.26) 

This operation looks naively like a resummation of the damping rate 7 in the propagator 
of the hard particles, but in reality it is the combined effect of a self-energy and a vertex 
resummation (recall the discussion after eq. (|2.9| ), and also the diagrams in Figs, ^.a and 
b); the resummation of the self-energy alone would have given an attenuation factor 27. 
The corresponding iterative solution reads then: 



KV - ' — 


f ■D + 7' 


w^(^) = 


7 / v-E 


f-D + 7\f-D + 7 


Kr(^) = 


7 /h/(^-i)\ 



n^ \-- /' (3.27) 

t> ■ iJ + 7 \ ' 

and suggests that, for very soft colour mean fields, the damping rate 7 ~ g^TXwiXj g) may 
act as an effective IR cutoff. That is, for gradients Dj. < g^T ^ we can even neglect the 
drift term v ■ D qjs, compared to 7 (at least, to leading logarithmic accuracy; cf. Sec. 3.4 
below), in which case the expansion ( ^.27| ) may be resummed into an exact solution (cf. 
Sec. 4.2). 

Diagramatically, the zeroth order solution W'^^^ in eq. ( p.27| ) corresponds to Fig. ^.b, 
i.e., to the sum of all the ladder diagrams in Fig. 0. (Incidentally, this is also equivalent 



to the relaxation time approximation, eq. (|1.6| ).) This is not a gauge-invariant subset of 



diagrams, and, indeed, it is quite obvious that the expansion eq. (|3.27| ) violates gauge 
symmetry at any finite order (since it treats the two pieces of the collision term on a 
different footing); see also the discussion at the end of Sec. 4.1. 

3.4 The leading-logarithmic approximation 

The previous applications of the Boltzmann equation (|2.15|) [0, |], ^ ^ have been mostly 
limited to the leading-logarithmic approximation (LLA) that we shall describe now. 

Recall first that, for colour excitations at the scale g^T ^ the collision term in 
eq. ( |2.15|) is known, strictly speaking, only to logarithmic accuracy, that is, up to correc- 
tions of 0(1/ ln(l/(7)) . This limitation has two sources: i) the IR problem of the damping 
rate |1^, |20|, ^, ^ , and ii) the gradient expansion in the presence of long range interac- 
tions m (i.e., the expansion in powers of S^^/g; cf. Sec. 3.2). Because of that, it has been 
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previously argued that the Boltzmann equation should be further simplified, for consis- 
tency, so as to preserve only the terms which are enhanced by a logarithm. Specifically, 
this involves two approximations: 



a) In eq. (|2.17|) for the collision integral $(v ■ v') one has retained only the singular 



piece of the magnetic propagator, namely p2| : 

rA(go«g)r^ ,^, \ ,,,, ^,^0 A^. (3.28) 



This allows one to isolate the IR singular piece of eq. (|2.17|) , which reads 



2 fv ■ v')^ 1 

<|.(v • v') ^ <|.o(v . v') = ^^ , ^ ^ In - , (3.29) 

TT^ml ./l _ (v • v')2 9 



where the logarithm ln(l/5f) in the r.h.s. has been generated via the following integral: 

-^^ = ln^^ln'-. (3.30) 

fi q /i g 

In this equation, the upper cutoff rriD is given by the screening effects at the scale gT (as 
included in *Vt, eq. ( p. 2^ )), while the IR cutoff /x is either the non-perturbative "magnetic 
p!7| (in which case /i ~ g'^T), or — in the framework of the effective theory for 



mass 



ultrasoft fields [0 — , the intermediate scale fi ^ g'^T ln{l/g) separating ultrasoft from soft 
momenta. In both cases, the estimate ( |3.30D holds to leading-log accuracy. By inserting 



the approximation ( p.29| ) into eq. ( p.l8[ ), we get the damping rate to the same accuracy 
(a = g^ATT) : 

7 ~ 70 = aA^Tln-. (3.31) 

9 



In fact, the expression of 7 obtained by evaluating exactly the integrals in eqs. (|2.17| ) and 



( |2.18|) with a sharp IR momentum cutoff fi (see Appendix B in the last paper of Ref. |]22| 



is: 

7 = «Arrin — , (3.32) 

/^ 

up to corrections of orderQ fi/m^,. 

b) The covariant gradient operator, or drift term, v ■ D^ ~ g^T in the l.h.s. of 
eq. ( |2.15| ) has been neglected next to the collision term oc 7 in the r.h.s. 

After these simplifications, eq. ( |2.15| ) reduces to: 



vE'^(x) = 7o{w^'^(x,v) - {W%x,y))o}, (3.33) 



''Actually, numerical studies of eqs. ( 2.17 ) and ( ^.18|) show that the error is even smaller, of order 
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where the subscript refers to the LLA, cf. eqs. ( p.29D and ( |3.31| ) (e.g., the angular 



average {W'^{x,y))q is given by eq. ( p.20|) with $ -^ $o)- 



Eq. ( |3.33| ) can be easily solved by iterations, as in eq. ( |3.27| ): the first iteration 
yields W^'^^ = v • E/70, and all the higher order iterations vanish (W^^^ = for A^ > 1) 
since W^^'^ is an odd function of v, while $o(v ■ v') is even (cf. eq. ( p.29| )). Thus, 



7o 

which, as already mentioned in the Introduction, is formally equivalent to the relaxation 
time approximation (|1.6|) , and generates a colour current j" = o"oE" , with the colour 
conductivity in the LLA ctq = m'jy/3'yo 0, 01- 

However, the approximation (b) above is insufficient for several reasons: 
i) It is incorrect in the electric sector, where it fails to provide Debye screening |^]. Indeed, 



eq. ( |3.34|) yields j° = 0, to be contrasted with the correct result ( |3.23| ): j° = —mjjAQ^x.). 



This is so since, for static fields, iy"(x, v) = —A'^{'x) is an exact solution of eq. (|2.15|) , 
for which the collision term vanishes (cf. eq. (|3.22| )); in this case, it is not legitimate to 
neglect the drift term. 

a) In some specific kinematical situations (essentially, for fields which are arbitrarily weak 
and slowly varying) , the linearized Boltzmann equation can be solved to a higher accuracy 
than in the LLA, leading to a formula for the transverse colour conductivity valid beyond 
the LLA. This will be explained in Sec. 4.2 below. 

in) In order to study the non-local structure of the ultrasoft amplitudes, one has to retain 
the drift term in the Boltzmann equation. Then, both pieces of the collision term (local or 
non-local in v, cf. eq. (|2.15|) ) play a role, as required by gauge symmetry (cf. eq. (|3.4| )). 
This will be further explained on the example of the polarization tensor, in the next 
section. 



4 The polarization tensor 

In this section we shall discuss in more detail the polarization tensor n^,y(P) for ultrasoft 
(P ^ g'^T) fields, as determined by the solution to the Boltzmann equation. The polariza- 
tion tensor typifies the non-local structure of all the ultrasoft amplitudes: as in the HTL 
approximation, all the n-point vertices with n > 3 follow from it as a consequence of the 
non-Abelian gauge symmetry (these vertices are generated by the covariant derivative in 
the l.h.s. of the Boltzmann equation (|2.15|) ). 
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4.1 Tensorial structure and iterative solutions 



According to eq. ( p.6| ), in order to construct tlie polarization tensor it is enough to consider 
a linearized version of the Boltzmann equation ( p.l9|) , namely: 

{v■^^)W%x,^r) = v■E%x)--f{W''{x,^r) - {W%x,^r))}, (4.1) 

where El{x) = 9M° — d^A]^ denotes only the "Abelian" piece of the electric mean field. 
It is then useful to go to the momentum representation, where eq. (|4.1| ) becomes 



{vP)W{P,v) = ivE{P) - t-f[w{P,v) - {W{P,v))}, (4.2) 

with P^ = {uj, p), E^{P) = i{ujA'^{P)—p'^A^{P)), and the colour indices have been omitted 
since trivial: the linearized equations ([4.1|) or ( ^4.2| ) are indeed diagonal in colour. 

Even though linear, eq. (|4.2| ) is still difficult to solve in general, since the term 
7(iy"(x, v)) is non-local in v (cf. eq. ( |2.20| )). Below, we shall consider an iterative 
solution, following the procedure explained at the end of Sec. 3.3. But before doing that, 
we shall derive from eq. (^4.21) some general properties of Il^^{P). 

First, the solution W{P, v) can be written in the form: 

W{P,v) = iW'{P,v)E\P), (4.3) 

with the new functions iy*(P, v) satisfying: 

{vP)W'{P,y) = v' - i-f[W'{P,Y) - (W^*(P,v))}. (4.4) 

The corresponding colour current can then be written as: 

jf'iP) = ml [^vf'W{P,v) = (y^\P)E\P), (4.5) 

J An 

with the following conductivity tensor: 

a^^'iP) = iml f^v^'W\P,v). (4.6) 

J An 

The polarization tensor is then defined by j^{P) = Il.^^{P)Ay{P), which implies: 

n^O(P) = -ip>a>'^{P), n'^*(P) = -iua^^^P). (4.7) 

By using eq. ( [4.4|) , one can verify that the resulting polarization tensor is transverse. 



Pfji^'^ = 0, and symmetric, U^'^ = U'^^. These properties, however, are not manifest on 
eqs. {^ and (|3). 
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The transversality property reflects the conservation of the (hnearized) current 
{P^j^ = 0) and has been already proven in a more general context in Sec. 3.1 (re- 
call eq. ( |3.7| )). Here, it immediately follows from eq. ( [1.41) , which implies (compare to 
eq. i3)): 



I "^ [v . P)W\P,w) = 0, (4.8) 

J 47r 



so that P^cr^* = 0, and hence P^II'^'' = 0. For what follows, it is useful to decompose 
iy*(P, v) into its longitudinal and transverse components with respect to p, by writing 
(with p* = pYp and p = |p|) : 

W\P, v) = Wl{P, ^)f + W4(P, v), Wl = V-^. P • Wy = 0, (4.9) 



and note that eq. ( |4.8| ) entails a constraint on the longitudinal component alone: 

|^(^_v.p)iy^(P,v) = 0. (4.10) 

By using this constraint, together with the properties of the angular integration, we shall 
now verify the symmetry property 11'^'^ = 11'^'^. Eqs. ( [1.6D and ( |4.7| ) imply, e.g., 

n°XP) = ujml j '^W^P.w), ff°(P) = pmlj'^v'WLiP.^). (4.11) 

These two expressions are indeed identical because: 

ffO(P) = fmlj'^{wv)WL = fuomlJ^WL = n°^(P). (4.12) 

In writing the first equality above, we have used the fact that p is the only remaining 
vector after performing the integral over v; then, eq. ( [4.1(j| ) has been used to obtain the 
second equality. It can be similarly shown that II*-^ = II-'*. 

The above properties fix the tensor structure of II^'^: as in the HTL approximation, 
n^*^ is determined by two independent scalar functions IIl{uj,p) and IIt{uj,p), which we 
choose as: 

n^(o;,p) = -n°0(P) = -pm^|^iy^(P,v), 

Uriu,p) ^ i(5,^._p,p^.)ff:'(P) = 1^^^|^vWt(P,v). (4.13) 

In terms of these functions, the components of 11'^'' read: 

n°°(P) = -n4^,p), U'\p) = -^Ul{u;,p), (4.14) 

rfi(P) = (S'l - p^fP)UT{u;,p) - ffP ^ nduj.p) . 
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As p — ^ 0, there is no privileged direction, and, since IIt{uj,p = 0) is non-zero (e.g., 
IIt{uj,p = 0) = mjj/S = Upi in the HTL approximation [jl^), the above expression for IT^ 
requires IIl{uj,p -^ 0) to vanish in the following way: 

2 

Ul{u,p^O) ^ -^nr(u;,p = 0). (4.15) 

The longitudinal and transverse components of the conductivity tensor will be also needed 
later. Writing j* = p^JL + Jt and E^ = p^Ei + E\, (with El = i{ujAL — pAo) and 
E'ij, = iuAlj,), and defining a^ and o"t such that Jl = ctlEl and j^ = (JtE^, we get: 

(Tl{uj,p) = -i — nLiuj,p), (Tt{uj,p) = -Ut{uj,p). (4.16) 

P^ 00 



We now turn to a discussion of the iterative solution to eq. ([4.4|) , following the 
considerations at the end of Sec. 3.3. If we treat the whole collision term as a perturbation, 
then the first two iterations read (cf. eq. ( |3.25| )): 

W^'^ = ^^ W^'^ = -, ^_ 1^^ _ /^^\1 (4 17) 

vP' ' vPXvP \vP/}- ^ ^ 

The first term above yields the well known HTL approximation for the (retarded) polar- 



ization tensor ||16| , [17[ , namely 



ni°j(c., P) = m^ -6X + - / ^ ,. !:':,.. \ ■ (4.18) 



An uj — v ■ p + ie 

The second term in eq. ( |4.17| ) gives then a correction to the HTL result which can be 
written in the following form: 

= -tUTTiL) ^ -— -—^{\-v) -<^ T^h (4.19) 



An J 4:71 V ■ P \v ■ P v' ■ P 



where in the second line we have used the definition (|2.20| ) of the angular averaging 
together with eq. ( |2.18| ) for 7. Higher order iterations HJ^^^ can be written down similarly. 



in a straightforward way. 



With the leading logarithmic approximation ( 3.29| ) for $(v • v'), eq. ( [4.19| ) coincides 



with the expression recently obtained by Bodeker in Ref. |3^ by diagrammatic calcula- 
tions. From the discussion in Sec. 3.2, one easily associates the first term within the 
braces in eq. ( [4. 19]) with the two diagrams in Figs. |^.a and b, and the second term, which 



is non-local in v, to the diagram with two hard loops in Fig. ^.c (the two unit vectors v 
and v' correspond to the velocities of the hard particles running around these two loops). 
The higher-order iterations H|j^^ with N > 2 would similarly correspond to the diagrams 
illustrated in Fig. |1^. 
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For P ~ g'^T , however, the contribution in eq. ( ^.19| ) is actually of the same order 
in g as the HTL ( 4.18| ), and even dominates over the latter by a logarithm hi{l/ g). This 



reflects the fact, already emphasized in Sec. 3.3, that the iterative expansion is generally 
not appropriate for the problem at hand, and it makes a priori no sense to try and evaluate 
11 ^jj from just a finite number of terms in this expansion. In the next subsection, we shall 
rather construct exact solutions to the equation (|4.4| ) in specific kinematical limits. 



Consider finally the second iterative solution, as described in eq. ( p.27| ). In Sec. 3.4, 
this proved to be useful in obtaining the leading-logarithmic estimate in eq. ( p.34|) . In 



general, however, this expansion must be used with caution since, as advertised at the end 
of Sec. 3.3, it violates gauge symmetry at every finite order. For instance, if we restrict 



ourself to the zeroth order iteration, we obtain (cf. eqs. (4^) and ( 3.27 )) 



lyf ) = -^—- , (4.20) 

which then leads to a polarization tensor which is neither symmetric, nor transverse. For 
instance, the conductivity tensor a*-' built out of ( [4.20| ) reads: 

.gV.p)=„n^/fl '•'"'^ , (4,21) 

which is not transverse in the static limit: p^ai^ [uj = 0, p) 7^ 0. 

4.2 Colour conductivities 

We now study the behaviour of the polarization tensor n^j,(c(j,p) at very small energy 
and momentum, cj, p <^ 7. This is interesting since the gradient expansion, which is only 
marginally justified for inhomogeneities at the scale g'^T , becomes more and more accurate 
as the inhomogeneity becomes softer and softer. (Strictly speaking, colour inhomogeneities 
cannot be unambiguously defined at extremely soft scales p <^ g'^T, for the reasons 
explained in the Introduction. The forthcoming discussion is nevertheless interesting since 
it applies also for momenta p ~ g'^T, at least within an expansion in powers of 1/ ln(l/(yf).) 
By solving exactly the (linearized) Boltzmann equations (|4.1| ) or ([4.4|) in this kinematical 
limit, we shall recover the previous result about Debye screening (cf. eq. ( |3.23| )), and 
compute the longitudinal and transverse colour conductivities defined in eq. ( [4.16|) . By 
"exact solutions" we mean here solutions which are obtained without using the leading 



logarithmic approximation ( |3.29| ) for $(v ■ v'), and which are known up to corrections of 
0{p /■-)). Of course, these solutions account only for the contributions of the hard and soft 
modes, with g > /i, to the corresponding conductivities. But they are still interesting as 
they allow for the matching with the corresponding contributions of the ultrasoft modes, 
to be computed non-perturbatively (see eq. (|4.36| ) below, and the discussion after it). 
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We shall study the two following situations: (i) u <^ p <^ 'y (this includes the static 
case L<j = as a particular limit), and (ii) p <^ u ^ •y. Since, generally, the electric and 
magnetic sectors behave differently in these limits, it is useful to project the Boltzmann 
equation ( [4 .41 ) for iy*(P, v) onto longitudinal and tranverse components (cf. eq. ( |4.9|) ): 



{vP)WL{P,^r) = vp - 27{iyz.(P,v) - (iyL(P,v))} 
{vP)W^{P,y) = 4 - i7{w4(P,v) - (H4(P,v))}. (4.22) 

{{) Consider first the static limit cu — > 0, where the longitudinal sector should provide 
Debye screening, as shown in Sec. 3.3. And, indeed, for u = 0, the above equation for 
Wl reduces to: 

{y-p){1+pWl) = ^l{WL{p,^r) - (W^l(p,v))}, (4.23) 

with the obvious solution^] Wl{uj = 0,j9, v) = —1/p. Note that the collision term in 



eq. ( |4.23|) vanish identically for this solution, which is therefore independent of 7 (and thus 
the same as in the HTL approximation). Since, moreover, W{u! = 0) = pWiiyj = 0)Aq{p) 
for static fields (cf. eq. (^73|)), this solution is clearly equivalent to W{u! = 0,p,v) = 
—Ao{p), as expected from the discussion in Sec. 3.3. When inserted into eq. ( [4.13| ), it 
yields: 

Ul{uj = 0,p) = ml, (4.24) 



which, together with the first equation ( 4.16|) , gives the behaviour of the longitudinal 



conductivity a^ at small u (and for arbitrary p): 



ml 



aL{u^O,p) ^ -tu^. (4.25) 



The results ([4.24| ) and ( |4.25|) are the same as in the HTL approximation |jT6|: at small 



frequencies, the electric sector is not affected by the collision effects (see also Ref. |3T[| ). 

In the same limit, however, important modifications occur in the magnetic sector. 
Consider, indeed, eq. ( [4.22|) for Wt at tu = 0: 



(vp)W^ = -4 + ty{w^{p,v) - (Vr^(p,v))}. (4.26) 

In the absence of the collision term (that is, to zeroth order in the iteration ( [4.17] )), this 
equation would imply W^ = — f^/(v ■ p — ie) (where the ie stays for retarded boundary 

*This solution can be obtained also from the iteration ( [4.1?j ), when written for Wl and uj ~ 0: then, 
the zeroth order term reads W£ {ut = 0) = —1/p, which, being independent of v, makes all the higher 
order corrections to vanish: Wj^ {uj ~ 0) for iV > 1 (cf. the discussion prior to eq. ([3.22)). 
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conditions, as in eq.( |3.11| )), which would then generate the small frequency limit of the 



HTL magnetic polarization tensor [16, 17 



uP{uj < p) ^ -i- uml / -— v\5^j - PiPj)v^ 5(v ■ p) = -i--m%, (4.27) 

thus yielding a^p {u <^ p) ~ {ti /2){m\)/p). These expressions, which are formally singular 
as p — s> 0, are correct only as long as p ^ 7. For ultrasoft momenta p ^ 7, they are 
modified by the collision terms, as we discuss now. 

Note first that, unlike in the electric sector where Wl{uj = 0) = —1/p is independent 
of V, here Wf(p, v) is a non-trivial function of v already in the zeroth order iteration. 
For such a function, the collision term in the r.h.s. of eq. ( ^.26| ) cannot vanish; it is thus 



a quantity of order 7, with respect to which the drift term in the l.h.s. of ( |4.26| ) can be 
neglected in the longwavelength limit p <^ 7. In this limit, the equation for Wt reduces 
to: 

4 = i7{w4(v) - (W4(v))}. (4.28) 

This is a simple equation which can be solved exactly. Specifically, since w^ is the only 
vector left in the problem, we can write: Wf (v) = [C ji^^v\^ with some coefficient C. 
Then, 

(ly^(v)) = (CA7)(4) = <Cli^)vlr. (4.29) 

where we have denoted Kt>* = (t>*), so that (cf. eq. (|2.20|) ) : 



/ f^$(vvO(v vO 






Physically, k = (cos a), where a is the angle made by the velocities v and v' of the 
colliding particles, and the brackets denote averaging with respect to the scattering cross 
section, cf. eq. (p.20|) ; obviously, |k| < 1. Then, eq. ( [4.28|) fixes the coefficient C as 
C = 1/(1 — k). To conclude: 



v^ 



W^t(v) = ---^ = -z ^, (4.31) 

7 1 — K 7 — 



with (recall eq. (pASD): 



5 = 7«: = m^ ^ / ^ $(v • v') (V . v'). (4.32) 

With the leading logarithmic approximation ( p.29|) for $(v ■ v') in eq. (|4.32|) , the angular 



integral over v' vanishes by parity. Thus, 5 is a finite quantity of 0{g'^T), which is 
completely determined by the present approximation. Its explicit evaluation, however, 
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requires the full expression ( |2.17| ) for $(v ■ v'). We write S = aNTS, and obtain S by 
numerical integration of eq. ( [4.32|) . The result is 5 = —0.20305024... . 



With W^{v) from eq. ( 4.31| ), it is finally straightforward to estimate the transverse 



conductivity, or polarization tensor, in the present kinematical limit (cf. eqs. ( |4.13|) and 



aTiuj = 0,p^O) = -ml 



— V • Wr(v) 

477 



m 



D 



CO, 



pi 



3(7 — 6) 7 — 5 



P'7-5 



(4.33) 



where tUpi = m^/v^ is the plasma frequency, that is, the frequency of the longwavelength 
{p -^ 0) collective excitations [l^, [l^]. Strictly speaking, eqs. (|4.33|) hold only for very low 
momenta p <^ g'^T; indeed, in their derivation above, we have neglected the drift term in 
the l.h.s. of ( [4.26| ), but we have kept the term 6 ~ g'^T coming from the collision integral. 
In the limit where ln(l/f7) is large, one can further expand these expressions in powers of 
5/7 ~ l/ln(l/f7) and get, to linear order. 



ariuj 



0) 




(4.34) 



where the neglected terms are down by, at least, two inverse powers of ln(l/5f). Re- 
markably, it turns out that, within the same accuracy, eq. (|4.34|) holds also for momenta 
p ~ g'^T, which is a case of physical interest^]. Indeed, since p/7 ~ 1/ \n{l/g) as well, one 
can solve eq. ( [4.26|) by iterations, as a formal expansion in powers of In^ = l/\n(l/g). 
This yields (compare to eq. ( [4.31| )) 



wUp,^) 



7 



1 + ^ 

7 



vp 

t 

7 



+ 0(ln 



.-2 



(4.35) 



which shows that p ~ g^T and 6 enter on the same footing in W^. But when constructing 
the colour conductivity, as in eqs. ( 4.13 ) and ( 4.16| ), the term in eq. ( 4.35| ) which involves 
V ■ p vanishes after angular integration, so we are left with the same expression for o"t as 
above, eq. ( [4.341 ). By using eq. ( |3.32| ) for 7, together with the numerical estimate for 6 



given above, we can rewrite eq. ( [4.34| ) in the following form: 

aNT 



0". 



T 



^uj-^pr^ g^T) 



OJ, 



■pi 



In 






(5 + In 



.-2 



(4.36) 



This expression, which represents the contribution of the hard and soft modes (with 
momenta k ^ jj) to the colour conductivity at the scale g^T, turns out to be useful for 



the matching with the corresponding contribution of the ultrasoft modes ||39 |. 



•'We thank Larry YafFe for this remark. 
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Eq. ( ^.36| ) , which holds to leading and next-to- leading order in an expansion in pow- 
ers of ln~^, extends the solution to the Boltzmann equation for aT{^ ^ p ^ Q^T) beyond 
the LLA of Sec. 3.4. Note that, even to this order, the transverse colour conductivity 
remains local (i.e., independent of the momentum p), as in the LLA. But this is only 
true within the accuracy indicated in eq. (|4.36|) : the first corrections to this equation, of 
0(ln~^), involve p^/7^ and thus are non-local. 



The expressions in eq. ( [4.33| ) should be also compared to the HTL result in eq. ( |4.27| ): 



this shows that it is essentially the damping rate 7 which cuts off the divergence of IIt, 
or o"r, as p ^ 0. This is in line with Drude's picture of the electric conductivity, and it 
is interesting to pursue this comparison even further, so as to emphasize the difference 
between colour and electric conductivity (see also Ref. [Q for a discussion on this point). 
The colour fluctuation induced by a uniform, and transverse, colour mean field E^, as 
determined by the solution above (cf. eqs. ( [4. 311) and (|4.3| )) : 



V ■ Et 11 
Vr(p<7,v) ~ = Tco^v-Et, with T^i n / \ ^ (4-37) 

should be compared with the charge fluctuation induced by an electric field, in the relax- 
ation time approximation |jl2[ (below, ttem = e^/An, with e the electric charge): 



Weiip) = reKp)v-ET, With tM —^ — ~ —-1—-. (4.38) 

aemln(l/aem)T^ e4Tln(l/e) 

Besides the loss of two powers of the coupling constant in the denominator {tcoi ~ 1/fl'^, 
as compared to r^i ~ l/e"^), the colour relaxation time appears to be independent of the 
momentum of the hard particles (unlike r^i, which is proportional to p). Thus, the present 
approximation for colour transport is formally similar to the relaxation time approxima- 
tion for electric charge, but with a shorter, and momentum independent, relaxation timeQ. 
We emphasize, however, that in the case of colour this is a property of the exact solution 
of the corresponding Boltzmann equation, and not a consequence of the "relaxation time 
approximation" . 



(m) Another interesting limit where the (linearized) Boltzmann equation ( [4.4] ) can be 
solved exactly is the longwavelength, or zero momentum, limit, p -^ 0, at fixed frequency 
uj (with cij ^ 7, for the present approximations to apply). Once again, an exact solution 
can be found because, once p is neglected, the vector fluctuation iy*(v) is necessarily 
proportional to v\ Moreover, if the momentum p is strictly zero, one cannot distinguish 
between longitudinal and transverse polarizations, so that Vr*(to',p = 0, v) = C{uj)v^ must 
be the solution to (cf. eq. ([4.4|)) : 



ujW\uj,y) = v' - i-f{W'{uj,v) - {W'{uj,v))}. (4.39) 



'^We thank Henning Heiselberg for a clarifying discussion on this point. 

33 



This is similar to the previous eq. ( |4.28| ), so that the corresponding solution reads (compare 
to eq. (HI)): 



W\u;,p = 0,^r) = —— -, (4.40) 

a; + 2(7 - 0) 



which yields the following, isotropic, conductivity tensor: 

Lu + i{'y - 5) 



^2 

a'^{uj,p = 0) = t6'^ -^ — - = 6'^a{uj,p = 0). (4.41) 



In particular, as cu -C 7, we obtain 



^2 

a{p = 0, cj < 7) ^ -^ , (4.42) 

7 — 



which is formally the same result as in the static case (cf. eq. ( [4.33|) ) except that it applies 



now to both the longitudinal, and the transverse, conductivities: <7l{p = 0) = arij) = 
0) = aijp = 0). This should be compared with the previous results in eqs. ( [4.25| ) and 
( |4.33|) : unlike the longitudinal conductivity, the transverse one appears to be continuous 
in the double limit uj —>■ and p —;► (in the sense that the two limits give identical 
results) . 

5 Conclusions 

In this paper, we have used the Boltzmann equation describing the relaxation of colour 
fluctuations in order to generate a set of gauge-invariant amphtudes for the ultrasoft 
fields, i.e., the fields with momenta ^ g^T. These amplitudes determine the response of 
the hard quasiparticles to longwavelength (A > 1/ g^T) colour mean fields. They define 
an effective theory for the ultrasoft fields, resulting from integrating out the modes with 
momenta larger than g^T in perturbation theory. 

The strategy which has been used in this paper, namely the use of the kinetic theory 
in order to construct effective amplitudes for the soft fields, is reminiscent of our previous 
construction of the hard thermal loops from coUisionless kinetic equations |l[ . The new 
element here is the inclusion of the effects of the collisions, which is essential since these 
are leading order effects for the colour excitations at the scale g^T. This results in two 
important differences with respect to the previous analysis of the HTL's: a) At a techni- 
cal level, kinetic theory appears to be the only workable approach toward the systematic 
construction of the ultrasoft amplitudes. Indeed, unlike the HTL's, which are one-loop 
amplitudes and have been originally obtained in a diagrammatic approach |jl8|, |T^, the 



ultrasoft amplitudes correspond to an infinite series of Feynman graphs which are conve- 
niently resummed, to the order of interest, by the solution of the Boltzmann equation. 
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b) As for their physical content, the main new ingredient in the ultrasoft aniphtudes is 
the effect of dissipation as a consequence of colhsions. In the leading logarithmic ap- 
proximation, and also in the next-to-leading order (NLO) approximation as defined in 
eq. ( |4.36|) , the dissipation is simply encoded in a local colour conductivity. In general, 
this is described by the (non-local) imaginary parts of the ultrasoft amplitudes (as, e.g., 
in eq. (glOD). 

By studying the Boltzmann equation, we have been able to obtain a few exact 
results about the ultrasoft amplitudes, in particular, the Ward identities they satisfy 
(cf. eq. (^3))' ^^6 static limit of the induced current (cf. eq. ( p.23| )), and the colour 
conductivity beyond the leading logarithmic approximation (cf. eqs. ( [4.33| ), ( [4.36| ) and 



( 4.41 )). More generally, by using formal solutions obtained by iterations, we have studied 
the non-local structure of the ultrasoft amplitudes (cf. eq. ( ^4.19| )), and established their 
diagrammatic interpretation (cf. Sec. 3.2). 

As already emphasized, the dynamics of the ultrasoft colour fields described by 
the Boltzmann equation is dissipative: any initial colour excitation will die away after 
a typical time Tcoi ~ l/{g'^Tln{l/g)). (This should be contrasted with the coUisionless 
dynamics in the HTL approximation, which is conservative [|^, and even Hamiltonian 
fill , pTJ.) This dissipative description is appropriate to study the relaxation of given 
initial off-equilibrium deviations in the average colour density, as, e.g., in the calculation 
of the colour conductivity, in Sec. 3.2. 

For many other applications — for instance, in studies of the baryon number viola- 
tion in the hot electroweak plasma — , one is interested in colour excitations at the scale 
g'^T which are generated by thermal fluctuations in the plasma. The most convenient 
strategy to deal with such non-perturbative fluctuations is to treat them as classical fields 
at finite temperature, which can then be simulated on a lattice [EB, B^, 53, |2^, 0, [2^]. In 



such a framework, one has to be able to also generate the correct thermal correlations 
of the ultrasoft fields, at least in the classical approximation. In practice, and especially 
for numerical simulations, it is convenient to use a Langevin description of the fluctua- 
tions, that is, to simulate the thermal correlations with an appropriate "noise" term. The 
"noise" is a random source with zero expectation value but non-trivial correlators which 
are chosen so as to induce, via the equations of motion, the proper thermal correlations 
of the ultrasoft fields. 

For the effective theory at the scale g'^T, the appropriate noise term has been con- 
structed by Bodeker [0. This term does not appear naturally in the derivation of the 
Boltzmann equation from quantum field theory [Q , where one focuses on the distribution 
function of the hard particles, rather than on the dynamics of the soft fields. Still, to the 
order of interest, the structure of the noise can be reconstructed from the known structure 



35 



of the collision term, by using the fluctuation dissipation theorem. This construction will 



be presented somewhere else [42 
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